First-order superfluid-Mott-insulator transition for quantum optical switching 
in cavity QED arrays with two cavity modes 
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We theoretically investigated the ground states of coupled arrays of cavity quantum electrodynam- 
ical (cavity QED) systems in presence of two photon modes. Within the Gutzwiller-type variational 
approach, we found the first-order quantum phase transition between Mott insulating and super- 
fluid phases as well as the conventional second-order one. The first-order phase transition was found 
only for specific types of emitter models, and its physical origin is clarified based on the analytic 
arguments which are allowed in the perturbative and semiclassical limits. The first-order transition 
of the correlated photons is accompanied with discontinuous change in the emitter states, not only 
with the appearance of inter-cavity coherence in the superfiuid phase. Our result gives a clear insight 
on the condition for the first-order transition to occur, and can give a strategy for future design of 
quantum optical switching devices with cavity QED arrays. 

PACS numbers; 42.50.-p, 42.50.Ct, 73.43.Nq, 71.36.-l-c 



I. INTRODUCTION 

Recently, parametric controls of quantum optics 
through manipulating states of matter in quantum op- 
tical systems have become key technologies to develop 
a new quantum optical device. In the studies, a cavity 
quantum electrodynamical (cavity QED) system, where 
photons interact with one or N emitters in one cavity, 
has been known to exhibit striking physics such as the 
strong nonlinearity, open-dissipative nature, and phase 
transitions In this context, coupled cavity QED 

arrays have also been focused so much in recent years, due 
to the new possibility toward a manipulation of the quan- 
tum optics with an efficient use of the many-bod y fe ature 
of photons after pioneering papers appeared jla - [l5l | . To- 
ward the realization, the effective coupling between cav- 
ity arrays was recently obtained with photonic crystal mi- 
crocavities [l6j , although its strong light-matter coupling 
regime has not been reached yet. The simplest case of 
the cavity QED arrays is described by Jaynes-Cummings 
Hubbard (JCH) model, where the photons can hop to the 
neighboring cavities while they suffers repulsive interac- 
tion with each other through the light-matter interaction. 

Coupled cavity QED arrays have been shown to ex- 
hibit the superfiuid (SF)-Mott-insulator (MI) transition 
of photons |12| - [l4j . being similar to Bose- Hubbard (BH) 
model which has long been studied for Josephson junc- 
tion arrays and neutral Bose gases in optical lattices [TtI - 
[T9| . Therefore, they have been considered as a new can- 
didate of quantum simulator of many-body physics in 
solids. As well as the scientific interests, the quantum 
phase transition of photons is important also on its ap- 
plication purpose for a new source of quantum-correlated 



photons. 

Recently, numbers of papers in this issue focus on fur- 
ther aspects on the pure physics e.g. the quantum fluctu- 
ations, ultrastrong coupling physics, and nonequilibrium 
physics (20l - l29j . On the other hand, the technological 
applications of cavity QED arrays have not been focused 
so much. For applications, the quantum phase transition 
of the coupled cavity QED arrays seems to be useful for 
optical switching [30[ and sensing devices with high sen- 
sitivity. However, in order to obtain a high performance 
of the switching, the phase transition of the first order 
will be more feasible than the conventional second-order 
one for the simple JCH model. Therefore, important 
questions to clarify for applications are now, (I) does the 
first-order transitions exist in cavity QED arrays, and if 
yes, (II) what is the condition for the first-order transi- 
tions to arise? In this regard, it is known in BH models 
(with a close relation to cavity QED arrays [3l|) that 
the SF-MI transition can be of the first order if the mul- 
tiple components of bosons are present [s^ - lssf . It was 
also shown in spin-1 bosons that the first-order transition 
is only found for antiferromagnetic spin-spin interaction. 
Therefore, both the multiple components of bosons and 
the types of boson-boson interactions will be the key in- 
gredients for the first-order transition to occur. 

In this paper, we study the SF-MI transition in a va- 
riety of coupled cavity QED arrays in presence of two 
cavity modes (Fig.[T]), as a simplest case of multiple com- 
ponent bosons. We find the first-order SF-MI transitions 
can occur also in cavity QED arrays and can be applica- 
ble to optical switching devices. The Hamiltonian used 
here is generally given by 



i— site z,j m—A.B 



(1) 
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where the first and second terms are the Hamiltonian 
within each cavity i and hopping of photons of modes 
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FIG. 1: Coupled cavity QED array with two cavity modes. 
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m (= A, B) between neighboring cavities i and j. The 
types of the hght-niatter interactions are varied with a 
choice of fuu. By applying Gutzwiher-type variational 
approach |36| for the ground states of coupled cavity 
QED arrays, we show a certain range of the models ex- 
hibit a first-order phase transition. Therefore, we con- 
clude that the types of light-matter interactions are im- 
portant for the first-order transition to occur also in this 
coupled QED arrays, being consistent with the case of 
spin-1 BH models jsiH]. 

This paper is organized as follows. In Sec. |IT1 we show 
detailed results of the Gutzwiller-type variational calcu- 
lation for a A- type three- level configuration of emitters as 
an example of the model in Eq. ([1]), where the first-order 
transition is found. In Sec. IIII| we show results for four 
different types of emitter models where the first-order 
transition is found only for specific types. Conditions for 
the first-order transition to occur are discussed. Finally, 
we will summarize the results, conclusions, and future 
remarks in Sec. IIVI 



FIG. 2: Energy diagram for A-type three level emitters. 
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FIG. 3: Lowest energy levels i.e. the eigenvalues of Eq. ((2|, 
for A-type three level emitters for A/g = 1.0. 



II. TWO-MODE CAVITY QED ARRAYS WITH 
A-TYPE THREE LEVEL EMITTERS (Al) 

In this section, we focus on the quantum phase transi- 
tions in two-mode cavity QED arrays with A-type emit- 
ters, as a simple example showing the first-order phase 
transitions. 



A. Dressed states of A-type system (t = 0) 

Wc will first show the detailed results for a two-mode 
cavity QED array system with A-type three level emit- 
ters whose energy level diagram inside a cavity is shown 
in Fig. [5] We call here this model "Al" . For this config- 
uration, photons of a mode A couples with a level tran- 
sition between |0) and |1), and a mode B couples with 
a level transition between |1) and |2). The Hamiltonian 
ho^i in Eq. ([T]) is given by 

ho = ujx\l){l\ + ^\'2){2\+uJAa\aA + ujBalaB 

+ (gio|l)(0|aA+5i2|l)(2|aB + h.c.)-MMot, (2) 



where a a and as represent annihilation operators of 
cavity modes A and B. The site index i is omitted 
for simplicity. The total excitation number per cavity 
Aftot = |1) (ll+a^aA+fls^s is conserved with this Hamil- 
tonian ho and the mean excitation number is fixed by a 
chemical potential ^. Throughout this paper, we assume 
that /U measures the strength of the energy injection by a 
external pump bath. Here we also assume 510 = gi2 = 9 
and a resonance condition for cavity modes llia = i^x and 
lub = (jJx — A for simplitity. We note that this simplifi- 
cation does not change the discussions. Without loss of 
generality, we set A > in this paper. 

Without tunneling t = in Eq. ([1]), the energy lev- 
els of this A-type cavity QED system are given by the 
eigenvalues of Eq. ^ and shown in Fig. |31 In the fol- 
lowing, we use a vector \a, ha, tlb) with an emitter state 
a (= 1,2,3) and numbers of photon modes ua and ub 
(eigenvalues of fiA (= a^a^) and ris (= a^os) equal to 
0, 1, • • • ) that forms the complete basis set of this Hamil- 
tonian. The eigenstates of the Hamiltonian ho consist 
of dressed states and decoupled states. There are three 



types of dressed states given by 

\M,nA,nB) = 



1 + ub 



2 + UA + riB 



l + UA 



UA + riB 



\0,nA + 1,71b) 



-\2,nA,nB + 1), (3) 



and 



|±,nA,ns) 



±\l -\l,nA,nB) 



l + UA 



I + ub 



4: + 2nA + 2nB 
with the eigen energies 



4 + 2nA + 2nB 
2,nA,nB + 1), 



\0,nA + l,nB) 
(4) 



^^|M,n.4.nB) {ujx ~ ^^)inA+■nB + 1) - AUB, (5) 
E\±.nA.nB) = {i^x - fJ'){nA + nB + I) - AUB 

(6) 
(7) 



±5 v^iX+n^B^rs . 

The remaining decoupled states are 

|0,0,nB) and |2,7m0), 

with the eigen energies 

E\o,o,nB) = {ujx ~ fJ. ~ A)nB, 
E\2,nA,Q) = {^x - fJ.)nA + A, 



(8) 
(9) 



Figure [3] shows that the ground state changes from the 
vacuum state |vac) = |0,0,0) to dressed photon state 
\ — ,0,nB) {riB = 0, 1,---) with larger total excitation 
number Aftot = J^b + 1 i^ftot is the eigenvalue of A/jot) 
as the chemical potential /z is increased. For ux ~ M ~ 
A > 0, the ground state changes at level crossing points 
fi ^ fic with (^c - ujx)/g = -V2, -A/g - {VS - V2), 
-A/g - (\/4 - V3), • • • -1.414, -1.318,-1.268, • • • 
for A/g = 1.0). 



B. Ground state phase diagram {t 7^ 0) 

Now we will find the ground state of the cavity QED ar- 
rays in case of t 7^ when the hopping of cavity photons 
is effective. Using a Gutzwiller-type variational wave- 
functions for the coupled cavity system 

m - n 1*^)' (10) 

'i— site 

= ^ g{a,nA,nB)\a,nA,nB), (11) 

a,nA.nB 

where the ground states of a product state with site in- 
dependent |$i). We minimize the expectation value of 
the free energy ($|7?|$) to find the variational param- 
eter g(a^nA,nB)- This approach gives the same result 
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FIG. 4: Ground state phase diagram of two-mode coupled 
cavity QED arrays with A-type emitters obtained by the 
Gutzwiller-type variational wavefunction approach for A/g = 
1.0. The order of the phase transition are indicated at the 
phase boundaries, (bold, black) and (solid, blue/red), by 
"1st" and "2nd", respectively. All the second-order phase 
boundary is given by a perturbation theory in Eq. (|14[l . 
whereas the dashed line (blue) inside SF(B) phase is a part of 
the curve given by Eq. (|14|l plotted for a guide. Dotted line 
(black) in SF(B) phase shows a prediction of the first-order 
phase boundary between MIO and SF(B) phases within the 
semiclassical approximation in Eq. (|21|l . The dotted line is 
shown only for the low-^ region near the MIO phase bound- 
ary where Eq. (|21|) is meaningful. 



as the mean field approach |12| which is to approximate 
Eq. ([l} by a mean field Hamiltonian 



i 



(12) 



and minimizing its lowest eigen energy with respect to 
site- independent variational parameters ipA and tpB de- 
fined by ipA = {aA,i), "0-8 = i&B.i)- In Eq. z denotes 
a number of cavities of nearest neighbors. 

In Fig. m the ground state phase diagram is obtained 
numerically with a cutoff n„iax (= 15) to the photon num- 
ber basis so that < nA,nB < 'T-max- In the black area 
determined by wx— /i— A— < 0, the number of photons 
diverges for nmax — > +00 and we can no longer determine 
a stationary state with a fixed mean- number [3l] |. The 
phase diagram contains three phases: MI phases with 
ipA = V'-B = and quantized number states oi Mtot, and 
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TABLE I: Coefficients in the second-order perturbation energy to MI states for A-type tliree level emitters in Fig. [S] 



MIO 



unperturbed state Ca.i ~ 



Cb.i — 
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FIG. 5: Mean field energies are plotted as a function of a 
variational order parameter xpB ~ (o-b) for (a) {fi — (jjx)/g ~ 
— 1.55 and (a) (/i — LJx)/g ~ —1.35. We set = in this 
figure. 

two SF phases, SF(A) with i/'A 7^ and tpB = 0, and 
SF(B) with ipA = and tpB 7^ 0. Since the effective 
mean field Hamiltonian is equal to Eq. ([T]) with < = in 
MI phases, the ground states arc the vacuum states |vac) 
with A/tot = or dressed photon states | — ,0,71^) with 
A/tot = riB + I- As indicated at the phase boundaries in 
Fig. 21 the phase transition can be of the second or first 
order. 

For high (high excitation density), the phase tran- 
sition is always the second-order one between MI with 
Atot > 1 and SF phases (A or B). The second-order 
phase boundary curves between MI and SF phases are 
also given by a second-order perturbation theory with 
respect to the small zttp's in Eq. ([T^ . The energy cor- 
rection per cavity {SEi) to the unperturbed MI states 
{I = MIO, Mil, MI2 • • • ) has the following form: 

+ {zt - zH'^Cba) IV-bP + O ((teV)^) , (13) 

where the analytic expression of the coefficients Ca,i and 
Cb.i are given in Table HI 

If Ca,i > Cb,1: there are three possible cases: (i) 
zt < 1/Ca.i where the coefficients of iV'Ap and iV'sP 
in Eq. are all positive, (ii) 1/Ca,i < zt < 1/Cb,i 



where the coefficients are negative for iV'Ap and positive 
for IV'Bp, and (iii) 1/Cb,i < zt where the coefficients of 
\'4'a\'^ and are both negative. This indicates that 

the ground state is the MI state (V'a ^ i-'B = 0) for (i), 
and SF(A) state (i/'a > and ^b = 0) for (ii). As for 
(iii), the perturbative expression in Eq. (fT3)) around MI 
state does not give any reliable prediction. Therefore, the 
second-order phase transition from MI to SF(A) occurs 
at zt ~ XjCA.i for Ca,i > Cb,i- We also find, by applying 
a similar argument for Cb,i > Ca.i, that the second-order 
phase transition from MI to SF(B) occurs at zt = 1/Cb,i 
for Ca,i > Cb,i- As a result, the phase boundary of the 
second-order transition from MI to SF(A)/SF(B) is given 

by 

zt = min(l/CA,;,l/CB,/). (14) 

On the other hand, for low /i, the phase transition 
between the MI phase with A/tot = and SF(B) phase, 
and that between SF(A) and SF(B) phases are first-order 
transitions in Fig. S) The first-order transition occurs 
with a discontinuous .jump in the variational parameters 
i.e. a discontinuous change in the ground state. There- 
fore we have to deal with a nonzero value of tp in or- 
der to determine the phase boundary which cannot be 
obtained by the second-order perturbation theory from 
the MI phases. Instead, we evaluated numerically the 
mean-field energy profile in (i/ja, ■0B)-plane which shows 
the distinct two energy minima in Fig. [5] (a), and the 
position of the global minimum changes from tpB ~ ^ 
to tpB ~ 1.6 at a critical hopping parameter t = tc 
(—0.547 < logiQ (ztc/.g) < —0.53). If the same applies to 
the case of the second-order transition (Fig. [5] (b)), the 
energy profile shows only one minimum that smoothly 
changes with hopping parameter. 

C. Classical and Quantum optics, and 
Photoluminescence spectra 

In Fig. [6l we show the amplitude ij.'m [(a-c)], mean 
numbers (rim) [(d-f)], and second-order number correla- 
tion at zero time delay gm (0) = (al^al^amam) / (fim)'^ 
(which is not defined for a vacuum state) [(g-i)] for cav- 
ity modes m = A,B. They are plotted as a function 
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of hopping parameter for different chemical potentials: 
{fi — ujx)/g = —1-36 (left row), —1.41 (middle row), and 

— 1.50 (right row). Depending on the chemical potential, 
we found different types of phase transitions: A second- 
order transition from Mil to SF(B) for (/i — ujx)/g = 

— 1.36, a second-order transition from Mil to SF(A) 
and a first-order transition from SF(A) to SF(B) for 
{fx—LVx )/g = —1.41, and a first-order transition from MIO 
to SF(B) for (/X — uJx)/g — -1.50. The signatures of the 
first-order transition are the discontinuous jumps in these 
figures, which can be applicable to switching devices of 
the classical and quantum optics. For example, the lumi- 
nescence from mode B will show a strong antibunching 

(2) 

{g^ {^) = 0) at logjg(z</,g) = -1.05-0 and almost per- 
feet coherence (5^ (0) w 1) at \og-^^Q{zt / g) = —1.05 + 
in Fig. [6] (h). As being the typical first-order transition, 
metastable states are found near the first-order phase 
boundary. In Fig. [5] (b), (e), (h), (c), (f), and (i), the 
dotted curves extended over the first-order phase bound- 
ary show the plots for the metastable states i.e. the lo- 
cal minima of the mean field energy, and the end of the 
dotted curves indicate where the metastable states have 
disappeared and the local minima turned into local max- 
ima. 

The calculated photoluminescence (PL) spectra are 
shown in Fig. [71 Two figures, (a) and (b), are obtained for 
the parameters {^—u!x)/g = — 1.4075±(5 with S = 0.0025 
on cither side of the first-order phase boundary between 
SF(A) and SF(B) phases. With such small change in 
parameters, the large difference is obtained, and we see 
the first-order transition can be used also for a classical 
optical switching. The PL spectra are obtained using 
the linear response theory by assuming that the photon 
of mode m (= A, B) leaks out of the QED arrays very 
weakly [s^], and given by 



(|i-cox)/g=-l-36 
Mil SF(B) 



cx Re / {G\al{T)a,n\G)e-''^^-'^^dT 
7r|(^|a™|G)p 



{UJ + El - EgY + {irf 



(15) 



where \G) (|/)) and Eq (Ei) are the eigen state and eigen 
energy of the ground (excited) state for the mean field 
Hamiltonian in Eq. (jl2p given by the solution of the vari- 
ational problem of ip^s. The Heisenberg representation 



is used for aj„(r) 



exp[iH^^T]al, exp[ 



phenomenological parameter, 7^, is introduced here to 
account for a resolution of the detector. In Fig. [7] (a) ob- 
tained in the SF(A) phase, the PL from the modes A and 
B shows double Lorentzian peaks. They can be under- 
stood as spontaneous emission. This is consistent with 

(2) (2) 
the observation that g\ {0) and 5^ (0) are much less 

than 1 in SF(A) phase in Fig. [6](h). On the other hand, 
in Fig.[7](b) obtained in SF(B) phase, a strong main peak 
and the subpeak structure are newly found at the chem- 
ical potential {uj = fi) and the lower energy side (w < /i), 
respectively. The strong main peak corresponds to the 
coherent emission from the Bose condensate of mode B, 
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FIG. 6: Amplitudes, tpA and Vs, mean numbers, (fiA) 
and (ns), and second-order correlation at zero time delay, 

g'^A (t = 0) and gg\T — 0), are evaluated as a function of 
hopping parameters. The dotted curves extended over the 
first-order phase boundary show the values for the metastable 
states (local minima of the mean field energy), and the end of 
the dotted curves indicate where the metastable states have 
disappeared. 



and the subpeak structure correspond to the MoUow-like 
side peaks. In the Mollow triplet in quantum optics [38j . 
which are found for resonant fluorescence from an emit- 
ter driven by classical field in an out-of-equilibrium con- 
dition, the side peaks are found in both the higher and 
lower side of the main peak. On the other hand, in our 
case, the whole system of photons and emitters are in 
a ground state of the thermal equilibrium. Due to the 
difference, the subpeak structure should be found only in 
the lower side of the main peak (See Appendix IX)) . 



D. Origin of the first-order phase transition 

Here we discuss the physics why the first-order transi- 
tion appears in this system (especially for small /i region). 

In this two-mode cavity system coupling with A-type 
emitters, we find only the two supcrfluid phase SF(A) 
and SF(B), while we find no SF(A+B) phase with tpA ^ 
and TpB ^ (which can arise in different emitter models 
as shown in the next section). We can naively under- 
stand why there is no SF(A-I-B) phase in this system by 
the following semiclassical analysis (which can be reli- 
able in presence of sufficiently large amplitude of pho- 
tons IV'ml ^ 1 inside SF phases) in combination with the 
complementary results based on the perturbation the- 
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(a) SF(A) 



(H-a,x)/g = -1.410 (b) SF(B) 
log,o(z//g) = -1.05 \ ' 



(|J-tOxVg = -1.405 
log,o(z(/g) = -1.05 
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FIG. 7: Photoluminescence spectra of Al model {A/g = 1.0), 
SA{ijj) and Ss((xj) in Eq. (jlSp . are shown by blue and red 
lines respectively, for different parameters in SF(A) and SF(B) 
of Fig. m which are located on either side of the first-order 
phase boundary at \og-^^Q{zt / g) = —1.05: (a) SF(A) phase for 
{fi — ujx)/g = —1.410, and (b) SF(B) phase for {^ — ujx)/g = 
— 1.405. For the plots, the detector resolution is 7r/(7 = 0.1. 



ory given in Sec. IIIBI In the semiclassical approximation 
(putting ipm into a,„ in Eq. ([T]) and Eq. ([2])), the lowest 
eigen energy -Escmi of the system is 

Esemi = {UJX - ^J. - zt)\lpAf + (l^X - ^-l - Zt - A)\'iPb\^ 
+ I^bP + Scmit + (16) 

where the last term -Eomit is a small contribution from 
the average of total energy of the emitters: 



emit 



= {ujx - + AP2, 



(17) 



with Pa = |Q^)(a| being a projection to the emitter state 
a = 0, 1, 2. In the semiclassical limit ji/'Ap + IV'sP ^ 1; 
it is explicitly given by 



E. 



'omit 



2 |^^|2 + |V-bP^ 



(18) 



which is small compared to other terms in Eq. (1161) . Ne- 
glecting the small contribution, the saddle point of -Esomi 
is determined by 



(ux - M - zt)iJjA 
{ujx - y.- zt- A)-tpB 



(19) 
(20) 



The semiclassical equations have three solutions: (i) 
ipA = '0-8 = corresponding to MI, (ii) ipA = 
{g/2)/ (lux — fi— zt) and — corresponding to SF(A), 
and (iii) ipA = and ipB = ig/'^)/{^x — n — zt ~ A) 
corresponding to SF(B). Among the three, SF(B) of (iii) 
minimizes the energy in Eq. (|16p for A > 0. Therefore, 
there is no possibility for the SF(A-I-B) phase to appear 
in the semiclassical regime (large t and/or high ^). 

Here we should note that the semiclassical result is con- 
sistent with the perturbation theory presented in Scc. lII B] 



(small t) which gives no possibility of the phase transition 
from MI to SF(A-|-B) phase. However, contrary to the 
semiclassical equations, the ground state can be SF(A) 
instead of SF(B), even if A > for quantum regime with 
low ^ i.e. weak excitation regime. The results in the two 
limits (semiclassical and perturbative) suggest that, at 
low the ground state will change from MI to SF(A) 
and eventually to SF(B) as increasing the hopping pa- 
rameter t. In this case, the phase transition from SF(A) 
with {iPa > 0,tpB = 0) to SF(B) with (0^ = 0,V'b > 0) 
should be the first-order one. This is because, if SF(A) 
and SF(B) connect smoothly by the second-order phase 
transitions, other phases, MI or SF(A-I-B), must exist be- 
tween SF(A) and SF(B), which seems to be unphysical. 

Moreover, another type of the first-order phase transi- 
tion is possible, which is the transition from the MIO state 
(vacuum ground state) to semiclassical SF(B) state. The 
situatioir occurs when i?scmi('/'A = 0,-05 > 0) = £^mio = 
at the minimum of Eq. p^ . Restoring the last term 
^-omit in Eq. p6p . this condition is given by 



cjJx - + A 



0, 



(21) 



2 (jjx — /i — A — 

and plotted by a dotted line in Fig. U showing a good 
agreement with the exact numerical solution (thick solid 
line) in the low-/i limit. 

Physically, the above mathematical arguments can be 
understood as follows. Because the photon energy of 
mode B is less than that of the mode A (wb = wa — A), 
the dressed photon energy -Eph, namely a sum of the 
first three terms in Eq. ([T6)) . is reduced for large ub 
state if ua + ub is fixed. (This is also the reason why 
"■B ^ n-A = in all MI phase with Aftot ^ 1 a-nd why 
only SF(B) phase appears in semiclassical regime.) Re- 
minding that the modes A and B couple, via the dipole 
transition, to the emitter state \a = 0) and |a = 2), re- 
spectively, the large ub state has a probability P2 = {P2) 
larger than Pq = (Pq). However, in the case of P2 > -Poi 
total energy of emitters -Eomit becomes large, hence can- 
not be minimized. Therefore, there can be two competing 
states minimizing -Escmi = -£'ph + -E-omit , which correspond 
to the two distinct minima: (i) one minimizes Ep\^ but 
maximizes -Eemit — which is related to SF(B) phase; (ii) 
the other minimizes -Eomit but maximizes Eph — which 
is related to MIO and SF(A) phases. The competition 
between two states becomes effective only where the to- 
tal excitation number is small (i.e. /x is small) such that 
A/tot < 0(1). This is due to the fact that Eph and i?omit 
become comparable only for small fj,. 

Validity of this interpretation is clearly checked in 
Fig. [8] showing the probability of the emitter states Pa 
obtained for the ground state as a function of the hopping 
parameter. Being consistent with the above physical ar- 
guments, we find Pq > P2 in MIO and SF(A) phases — (i), 
whereas Pq < P2 in SF(B) phase — (ii). These two class 
of states does not connect to each other with a small per- 
turbation, and showing the first-order transition between 
(i) and (ii). We should also note here that the Mil, MI2, 
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FIG. 8: Probability of emitter states, Pa = {$|Pc|$) with 
Pa = are evaluated for the ground state |$) and plot- 

ted as a function of hopping parameter for different chemical 
potential: {ii — ujx)/ g = (a) —1.36, (b) —1.41, and (c) —1.50. 



and the MI states with higher A/tot smoothly connected 
with SF(B) state via the second-order phase transition. 
This can be understood by seeing P2 > Pq for the MI 
state | — ,0,713), since P2/P0 = 1 -I- > I for > 
from Eq. This means the higher MI states can be in 
the same class (ii) as SF(B) state. 

The origin of the first-order transition in our model 
is analogous to that of the first-order transition occurs 
in BH models with two or three component bosons [32I - 
[35} . In spin-1 BH model [l^l, the first-order transition 
occurred between two distinct MI and SF states: one 
is a MI state that minimizes the total spin-spin interac- 
tion energy but maximizes the kinetic (hopping) energy, 
and the other is a SF state that minimize the kinetic 
(hopping) energy but maximizes the total spin-spin in- 



teraction energy. The first-order transition is possible 
because, not only the intersite coherence but also the 
spin state (the internal degree of freedom) show a dis- 
continuous change at the phase transition. According to 
the discussion here, a presence of such competing effect 
showing distinctive multiple energy minima, due to an 
additional degree of freedom, seems to be a condition to 
find a first-order MI-SF transition in general model. 



III. 



TEST FOR OTHER EMITTER MODELS 



Here we discuss the ground states obtained for four 
different types of emitters (shown in Fig. As will be 
shown here, some models exhibit the first-order phase 
transition while some docs not. Combining the results 
with those shown in Sec. [Hi we refer to the condition for 
the first-order MI-SF transition to occur in cavity QED 
arrays. 



A. Four sample models of emitters 

Four types of the emitters inside a cavity considered 
here are shown and labelled by QDl, A2, QD2, and TLS 
in Fig. ini They are explained as follows. 

(QDl) A quantum dot four level system with a biexci- 
ton state |XX), two exciton states jX^^) and |Xj,), and 
ground state |G), which is coupling with two cavity 
modes A and B as shown in Fig. [5] (a). Dipole transi- 
tion is assumed to occur only in case of the resonance; 
a transition between the biexciton and exciton states is 
coupled to the mode B {u>b = — U), and that be- 
tween the exciton and ground states is coupled to the 
mode A [uja = i^x)- For this model, the Hamiltonian 
ho^i in Eq. ([1]) is given by 



ho = u}Aa\aA+ ujEa^'gCiB + {2LJX - ^)\XX){XX\ 
+UJX {\X^){X^\ + \Xy){Xy\) - ^iMtot 



+gaB\XX){{X,\ + {Xy\)+h.c 
+gaA{\X,) + \Xy)){G\+h.c., 



where 



Mot = 2\XX){XX\ + \X,){X,\ + \Xy){Xy\ 



~a\aA 



(22) 



(23) 



(A2) A-type three level emitters are coupling with two 
cavity modes as shown in Fig. [9] (c) , which are different 
from the one discussed in Sec. [H] (Fig. [5]). For this model, 
the Hamiltonian /io,i in Eq. ([IJ is given by replacing the 
light-matter coupling terms in Eq. ^ by 

g{aA + aB)|l)(0| + gi&A + aB)|l)(2| + h.c. (24) 

(QD2) A quantum dot four level system, which is cou- 
pling with two cavity modes A and B with different man- 
ners from (QDl), as shown in Fig.[n](c). Both the modes 
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| Xx>/ \|Xv> 
mode A 



||G>/ 



(c)QD2 



mode A |xx> 
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\ Xx)/ \|Xy> 

J|G>1 



(b)A2 



cox-A 



I \ y mode A 
mode B 



\o)l 



»''\ |2> 



(d) TLS 

li> 



cox-A 



mode A 



lo> 



mode B 



We find the first-order phase transition in QDl [Fig. [10] 
(a)], while the phase transition is always of the second 
order in A2 [Fig. M (b)], QD2 [Fig. M (c)] and TLS 
[Fig.[TU](d)]. It should be noted that for the latter cases, 
we found ^pA 7^ and 'ipB 7^ everywhere in the su- 
perfluid phase, which is quite different from the result 
for the A-type emitters in Sec. [H On the other hand, 
the first-order transition from SF(A) to SF(B) is found 
for QDl being similar to a model Al, the A-type emit- 
ters in Sec. |3l According to the result, we conclude that 
the first-order transition can be found when the different 
modes are coupling to the different transition processes. 

Here, let's see the origin of the different results between 
(Al, QDl) and (QD2, A2, TLS) by applying semiclassical 
arguments to the A2 model being similar to that was 
shown in the previous section. In the semiclassical limit 



I "05 1 ^ I7 the energy of the system is given by 



-V2.9|^a + V'b|+0(1), (28) 



FIG. 9: Models of three types of emitters considered in 
Sec, mil Here we set frequencies of the cavity modes A (= uja) 
and B (= ujb) as uja = oJx and ujb = cjx — A in the same 
way as in Sec.HTl 

A and B coupled to all the dipole transition irrespec- 
tive of the resonance condition. This model would be 
more suitable than the model (QDl) for the case where 
the biexciton binding energy A is of the same order of 
magnitude as g. For this model, the Hamiltonian ho_i in 
Eq. ([T]) is given by replacing the light-matter coupling 
terms (the third and fourth lines in Eq. ([H])) by 

gi&A + aB)|XX)((X,| + {Xy\) + h.c. 

+giaA + aB){\X,) + \Xy)){G\ + h.c. (25) 

(TLS) A two level atom, which is coupling with two 
cavity modes A and B through the dipole transition, as 
shown in Fig. ^ (d) . For this simple model, the Hamilto- 
nian /io,i in Eq. ([T]) is given by 

ho = UJAal^aA+UJBa^B^'B +^x\i){M- fjMtot 

+5(aA + aB)|l)(0|+h.c., (26) 

where 

Ntot = +a\a A + alaB. (27) 



B. Ground state phase diagram {t 7^ 0) 

In Fig.[Tni we show ground state phase diagrams which 
are obtained by applying the variational approach in 
Sec. H] to the four models, QDl, A2, QD2, and TLS. 



and the saddle point equations in the semiclassical limit 
are 

2{ujx - M - zt)ipA - V2g = 0, (29) 
2{ujx - f^- zt- A)tPB - V29 = 0, (30) 

where the relative phase between ipA and tJjb is locked to 
zero for the energy minimization. Therefore, in the semi- 
classical regime, only the SF(A-t-B) phase with nonzero 
TpA and nonzero ipB is possible, in contrast to the case of 
(Al, QDl). 

The same conclusion is drawn also in the perturbative 
regime where ^ 1- If the perturbation 

theory is applied to the three models (QD2, A2, TLS), 
one obtains the second-order energy correction to the MI 
states in a general form 

SEii^l;A,^B) = zt\4>A? + zt\-^B? ~ zH'^CaMa? 

-zH^CbMbI^ - zh^CABAi'Ars+rAi'B) (31) 

with Cab.i 7^ [which vanishes for Al and QDl as found 
in Eq. (|13p ]. This perturbative expression gives only the 
possibility that the second-order phase transition from 
MI to SF states occurs for two modes at the same tran- 
sition point given by 

1 — ztCA,i —ztCAB.i _ Q 
-ztCAB,i 1 - ztCB,l 

This indicates that all the Mott phase is surrounded by a 
SF(A+B) phase. Therefore, reminding the semiclassical 
result again, the SF(A-I-B) phase shows a smooth connec- 
tion between the perturbative and semiclassical regime, 
being in contrast to the models (Al, QDl). According 
to the considerations, it is reasonable to conclude that 
there is no first-order phase transition in the three mod- 
els (QD2, A2, TLS). 
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FIG. 10: Ground state phase diagrams for three difFerent models for emitters: (a) a quantum dot model (QDl) with A/g = 0.75 
and its enlarged view (right panel), (b) a A-type emitters (A2) with A/g = 1.0, (c) a quantum dot model (QD2) with A/p = 1.0, 
and (d) a two-level atom model (TLS) with A/g = 0.5. We set a cutoff for photon number so that nA,nB < rimux = 16 for 
(a), riA + nB < imax = 20 for (b), ua + ns < ^max = 15 for (c) and (d). The black region is where the number of photons 
diverges where uox — ^i — A — zt<0 



IV. CONCLUSIONS 

We investigated the possible quantum phase transition 
in coupled cavity QED arrays in presence of two cavity 
modes. Within the Gutzwiller approximation, which is 
equivalent to the mean field approximation, the ground 
state phase diagrams were obtained for types of emitter 
models with various light-matter interaction configura- 
tions, and we found the first-order phase transition be- 
tween Mott insulator (MI) phase and superfiuid phase (or 
between different superfiuid phases) for a certain range 
of the emitter models. It seems to be necessary, in ad- 
dition to the multiple components of photons, for the 
first-order phase transition to appear that the different 
photon modes should couple with different transitions 



between emitter states separately. (A finite difference 
of the transition energies A 7^ seems to be important 
as well.) In that case, the perturbative regime near the 
Mott insulating phase and the semiclassical regime can 
have different types of the ground states, suggesting the 
first-order phase transition can occur in their intermedi- 
ate regime. 

The origin of the first-order phase transition is anal- 
ogous to that of spin-1 BH model, where the first-order 
transition occurs between the states with large difference 
in the internal spin degree of freedom; MI states includes 
only the lowest spin state, meanwhile the SF states in- 
clude the high spin states. Therefore, the kinetic energy 
and spin-spin interaction energy gives the two compet- 
ing minima of the total energy in spin-1 BH model. In 
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case of the two-mode cavity QED arrays, the competition 
between the dressed photon energy and emitter energy 
give rise to the first-order phase transition. The argu- 
ment presented here on the origin of the first-order phase 
transition should be quite useful for future application of 
these system to possible quantum optical switching de- 
vices. 

As a future remark, the quantum phase transition of 
the first order, which is considered here in the ther- 
mal equilibrium condition, can occur also in the out- 
of-equilibrium condition. We know the thermodynamic 
phase transition between the different ground states oc- 
curs in order to minimize the total free energy of the 
system of interest. Is there any quantity that can charac- 
terize the phase transition even in the out-of-equilibrium 
condition? In the community of nonequilibrium statisti- 
cal physics, it has been conjectured that the entropy gen- 
eration rate of the system of interest should be minimized 
in the nonequilibrium stationary state, even though this 
claim has not been clarified yet [s^. Therefore, the en- 
tropy generation rate might be such a quantity as to char- 
acterize the nonequilibrium phase transition also in the 
coupled cavity QED arrays. The problem is interesting 
also in case of the phase transition of the second order 
in single-mode cavity QED arrays, since the relation be- 
tween thermodynamic and non-equilibrium phase transi- 
tions has been a recent hot topic. The related problem is 
also found in cxciton-polariton systems in a semiconduc- 
tor microcavity both in experiments [40| and theory [4l| - 
l43j . Theoretically, the cavity QED arrays considered 
here are very useful to investigate this problem, since 
the density matrix of the system of interest can be fully 
solved within the mean field approximation irrespective 
of whether the systems are in the thermal-equilibrium 
or out-of-equilibrium conditions [1^ . The application of 
the cavity QED arrays to that problem will be discussed 
elsewhere. 
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FIG. 11: Schematic explanation for Mollow-like PL spectra 
the non-equilibrium and thermal equilibrium situation. Here, 
(n -|- l)-th and n-th dressed states are focused for a situation 
where a two-level system is driven by a coherent field V ~ 
^/n > 1. 

In presence of a coherent field ^ ^/n 3> 1 , the energy 
level of the two-level atom suffers an optical Stark effect 
and forms a ladder of dressed states. Consider the (n+l)- 
th and n-th dressed states, as in Fig. [TTJ Such level 
diagram is often used to account for the Mollow triplet 
structure of the resonance fluorescence where an atom 
is driven by a classical field ip. The initial state, before 
the decxcitation and fiuorescence, is found both in the 
two branches of the dressed states with (n -I- 1) quanta 
in the driven and out-of-equilibrium system as shown in 
Fig. [TT] (a). This situation gives conventional Mollow 
triplet spectra (the lower panel). In our case, however, 
the QED system is driven while it is kept in the thermal 
equilibrium and occupies the ground state. In this case, 
the initial state is the unique ground state and found only 
in the lower branch of the dressed states (so-called lower 
polariton branch) of {n + 1) quanta as shown in Fig. [TT] 
(b). In this case, the resulting PL spectra should display 
the main peak and only low-energy side peak (the lower 
panel) . 



Appendix A: Mollow-like side peaks in 
non-equilibrium and thermal equilibrium situations 



Here, we demonstrate that, when a cavity QED sys- 
tem is in the thermal equilibrium ground state and large 
coherent field is present, photo luminescence exhibits 
Mollow-like side peaks only in a low-energy side of the 
main peak. This is in contrast to the non-equilibrium 
case of conventional Mollow triplet where an emitter is 
driven by a coherent field. To simplify the discussion, we 
consider a two-level atom as an emitter. 
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